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Flavor oscillations of neutrinos are analyzed in the framework of Brans-Dicke 
theory of gravity. We find a shift of quantum mechanical phase of neutrino pro- 
0\ \ portional to GnAtu 2 and depending on the parameter u>. Consequences on atmo- 

spheric, solar and astrophysical neutrinos are discussed. 
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1 Introduction 

Among all alternative theories of gravity, the Brans-Dicke (BD) theory |3J provides the 
most natural generalization of General Relativity It can be thought of as a minimal 
extension of Einstein theory in which Mach's principle and Dirac's large number hy- 
pothesis (see, for example, f2|) are properly accommodated by means of a nonminimal 
coupling between the geometry and a scalar field 0, the BD scalar. The scalar field rules 
dynamics together with geometry and, furthermore, induces a variation of the gravita- 
tional coupling with time and space through the relation G e ff = 1/0. The gravitational 
constant Gn is recovered in the limit — > constant. Some recent experiments || seem 
to confirm a variation of the Newton constant on astrophysical and cosmological sizes 
and time scale. 

The effective action describing the interaction of the scalar field nonminimally 
coupled with the geometry and the ordinary matter is given by M 
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where R is the scalar curvature, C m is the matter contribution in the total Lagrangian 
density. The constant uj is determined by observations and its value can be constrained 
by classical tests of General Relativity. The consequences of BD action fll.lj ) have been 
analyzed for the light deflection, the relativistic perihelion rotation of Mercury, and the 
time delay experiment, resulting in reasonable agreement with all available observations 
thus far provided uj > 500 ||. On the other hand, bounds on the anisotropy of the 
microwave background radiation give the upper limit u) < 30 ||. Einstein's theory is 
recovered for u — > oo. In this limit, the BD theory becomes indistinguishable from 
General Relativity in all its predictions. 

Understanding if the BD theory of gravity may be considered as the right general- 
ization of Einstein gravity and, as a consequence, how it affects physical phenomena is 
an important matter. In this paper we will face this issue by considering neutrino os- 
cillations, calculating, in particular, the contribution to the quantum mechanical phase 
mixing induced by the non-standard coupling between the geometry and the scalar field. 
As we will see, such a correction does depend on the value of the parameter uj. 

It is well known that the problem of neutrino oscillations is still open, and the research 
of new effects in which they could manifest is one of the main task of the today physics. 
For this reason, the quantum mechanical phase of neutrinos propagating in gravitational 
field (usually the Schwarzschild or Kerr field) has been recently discussed by several 
authors (see H-[n| and references therein), also in view of astrophysical consequences. 



More controversial is the debate concerning the red-shift of flavor oscillation clocks, 
given by a term proportional to 

G N Am 2 M r B 
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first derived by Ahluwalia and Burgard || in the framework of the weak gravitational 
field of a star, with mass M. Here Am 2 is the mass-squared difference, Am 2 = |m|— mf |, 
E the neutrino energy, ta and rs the points where neutrinos are created and detected, 
respectively. They also suggest that the oscillation phase (|1.2| ) might have a significant 
effect for supernova explosions due to the extremely large fluxes of neutrinos produced 
with different energies, corresponding to the flavor states. 

This result has been confirmed in the paper by Grossman and Lipkin ||, and it 



has been also derived by Konno and Kasai [ I I ] under the assumption that the radial 
momentum of neutrinos is constant along the trajectory of the neutrino, strengthening 
the correctness of the Ahluwalia-Burgard arguments. Nevertheless, assuming that the 
neutrino energy is constant along the trajectory, Konno and Kasai show that the term 



]2|) is cancelled out, recovering in such a way the result of Refs. |fLz|, |13} |. 
Without pretending to solve or face here this controversy, which goes beyond our aim, 
this paper is a straightforward extension of the calculations of Ref. |E| in the framework 
of BD theory. In [fP2 |, the neutrino oscillation formula in a gravitational field is based 



on the covariant form of the quantum phase that arises due to the assumed mixing of 

massive neutrino. The result (i.e. the cancellation of G^Am 2 term) is the same of Ref. 

T3| , but it is derived without invoking the assumption that underlying mass eigenstates 



are emitted at different time. We find that the scalar field in ( |1 . 1|) nonminimally coupled 
to the scalar curvature induces a red-shift of flavor oscillation clocks in the quantum 
dynamical phase, which is proportional to 

G N Am 2 1 r B 

= log — . (1.3) 
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It vanishes in the limit tu — > 0[|. Eq. (|1.3| ) can be seen, in some sense, as a further test, 
in addition to the standard ones above discussed, for establishing the validity (or not) of 
the BD theory. 

The layout of this paper is the following. In Sect. 2 we shortly recall the Schwarzschild- 
like solution coming from BD field equations, which describe the static and stationary 
gravitational field generated by a mass M, and the corresponding expressions in the 
weak field approximation (for details, see the paper |1[]). Sect. 3 is devoted to the calcu- 
lation of the quantum mechanical phase for propagating neutrinos in the BD geometry. 
Conclusions are drawn in Sect. 4. 



1 The extension of the paper M (or M) to the BD theory does not give appreciable correction to the 
quantum dynamical phase. In fact, the corrective factor is of the form 

3 + 2w (1.4) 
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which is ~ 1 for uj ~ 500 and u) < 30. 



2 Static Spherically Symmetric Field in BD Theory 

Variation of the action ( |1 . 1|) with respect to the tensor metric g^ v and the scalar field 
yields to the field equations PJ 
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for the geometric part, and 
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for the scalar field. □ is the usual d'Alembert operator in curved space-time and T^ v 
is the momentum-energy tensor of matter. The line element describing a static and 
isotropic geometry is expressed as 



ds z 
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where the functions a and f3 depend on the radial coordinate r. The general solution in 
the vacuum is given by 
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where the constants, appropriately chosen, are given by 
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In the weak field approximation, the components of the tensor metric, g^ u ~ r]^ v + h^, 



reduces to the form H 
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As discussed in Introduction, the weak-field solutions ( |2.8| )-( |2TTT| ) have been analyzed 
for gravitational red-shift, the deflection of light and perihelion of Mercury (to be more 
precise, the last one requires an approximation up to the second order in M/r). In 
the next Section we will investigate the phenomenological consequences of BD solutions 
(p.8|) ( pTT0| ) on neutrinos propagating in such a geometry. 



3 Neutrino Oscillations in BD Geometry 

The effects of gravitational fields on the quantum mechanical neutrino oscillation phases 
have been analyzed in the semi-classical approximation, in which the action of a particle is 



considered as a quantum phase |I5| . In calculating such effects induced by BD geometry, 
we will use the same approximation. 

A particle propagating in a gravitational field from a point A to a point B, changes 



its quantum mechanical phase according to the relation [[15 
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where p M = rag ^(dx" / ds) is the four-momentum of the particle and ds 2 = g^dx^dx 1 ' . 
Following Ref. [O, the quantum mechanical phase becomes 
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Inserting the momentum of the particle, coming from the shell-condition, g^p^p" = m 2 , 

p r = e p ~ a VE 2 - m 2 e 2a (3.3) 

into Eq. (|3.2j ), and using the fact that dt/dr = e ,3 ~ a , one gets a difference phase A$ 
given by 
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(3.4) 
By using the weak field approximation, Eqs. ( |2.^ )- ([2Tl0| ), one can separate out the 



gravitational contribution to the neutrino oscillation phase, so that Eq. fl3T4|) can be cast 
in the form 

A$ = A$ + A^ , (3.5) 

where (restoring the constants c and h) 
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which represents the standard phase of neutrino oscillations, and 
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In deriving Eqs. ( |3.(j| ) and ( ^.7| ), we have considered ultra-relativistic neutrinos, E » m, 
where E is interpreted as the energy at the infinite (see |12| for details). The integration 
has been performed along the light-ray trajectory where E is constant. 

It is convenient to rewrite the phases ( |3.6|) and ( |3.7| ) in the following way 

o Am 2 MeV rrj-n . 

A$ « 2.5 • lO 3 -^- — - 4—^ , (3.8) 
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where M & is the solar mass. Estimations of the difference phases ( |3.8|) and ( |3.9|) are 
carried out for solar, atmospheric and astrophysical neutrinos. To this end, we will 
introduce the ratio q defined as 

= A$ 1 M \og{r B /r A ) 

q A$ ~ ' 2 + uj M Q (r B -r A )/Km- l ' ' 

q does not depend on the squared-mass difference Am 2 and on the neutrino energy E. 
For solar neutrinos, we use the following values: M ~ M , r A ~ r Earth ~ 6.3 • 10 3 Km, 
and r^ ~ r^ + .D, where D ~ 1.5 ■ 10 8 Km is the Sun-Earth distance. Eq. ( 3.10| ) gives 
the result 

g~10' 8 -^, (3.11) 
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which is an irrelevant correction to the difference phase ( |3.8|) . Analogous conclusion holds 
for atmospheric neutrinos. 

Concerning the astrophysical neutrinos, the effect could be more relevant and could 
be measured by terrestrial experiments. In fact, setting r B = ar A , 1 < a < oo and using 
the typical values of neutron stars, M ~ 1.4M Q and radius r A ~ lOKm as in Ref. 0, we 
get 

2 + uj a — 1 

Till now, our analysis has been done for radially propagating neutrinos. In the case of 
motion transverse to the radial propagation and near to the detection point r A we have, 
following ||, 

= Am'cCM r B - r A _J_ W MeV M_ rj^-r, 
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Then, the ratio between the difference phases ( |3. 13| ) and (|3.8|) is 
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For the numerical constants corresponding to Sun and Earth, we have 

± 1.5 ■ IP" 5 ± 5 ■ 10- 10 
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Using the above values for a neutron star, Eq. ( |3.14j ) gives the result 

«r L ~^-- (3-16) 

As discussed in Introduction, experimental data imply that the parameter u can assume 
the value uj > 500. For the lower limit, one gets from Eqs. (|3.12|) and ( |3.16| ), 



g~4-10- 4 ^^, g ± ~4-10" 4 , (3.17) 
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giving a correction of the 0.01 percent. Values of uj < 30, coming from the anisotropy 
of microwave background radiation, allow to get corrections of few percents, as one can 
immediately derive from Eqs. (|3.12j) and ( |3.16| ). Such contributions to the quantum 



mechanical phase of neutrinos, are very significant and could be considered as a test for 
establishing the validity of BD theory. 

4 Conclusions 

In this paper, we have analyzed neutrino oscillation phenomena in the framework of 
BD theory. We have derived a correction to the standard difference phase of the order 
Gat Am 2 , which vanishes in the limit uj — > oo, when the BD theory reduces to General 
Relativity. 

Estimation of such a correction has been carried out assuming for the parameter uj the 
values uj ~ 500 and uj < 30. Such values may be relaxed considerably with the advances 
in technology associated with astronomical observations and astrophysical experiments, 
making our corrections as a mean to discern between BD theory and Einstein's theory, 
in addition to the ones discussed in the Introduction. 

Nevertheless, BD is a particular case of scalar tensor-theories where one assumes 
that matter acts as source of scalar field 0, which generates the curvature of space-time 
associated to the metric. The strength of the coupling between the scalar field and gravity 
is determined, in these theories, by the function u; (</>), which is constant in the BD theory. 
Besides, a self-interaction potential V(<p) can be also introduced, generalizing in such a 
way dynamics of the field. 



The dependence of the parameter u on <fi could have the property that, at the present 
epoch, and in weak field situations, the value of the scalar field (fi is such that to is very 
large, leading to theories almost identical to General Relativity today, but for past or 
future values of 0, as in strong field regimes as for neutron stars, u could take values that 
would lead to significant differences from General Relativity. In this sense, scalar-tensor 
theories are richer than BD theory and could play a relevant role in the neutrino oscillation 
physics (and in Pound-Rebka or COW experiments, as well as in atomic systems in 
linear superposition of different energy eigenstates). This because the variability of the 
parameter u implies that, in some epoch, its value could be very small, and, in such 
a way, a correction to the quantum mechanical phase of 10% can be obtained (in this 



particular case, the factor (|1.4j) reduces the Ahluwalia-Burgard result to 15%, instead of 



20% as derived in Ref. |J). In a forthcoming paper we will face these issues. 

The authors would like to thank the referees, and in particular D.V Ahluwalia, for 
the useful comments on the subject treated in this paper. 
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